Tetrahedral and Triangular Deformations of Z = N Nuclei 

in Mass Region A ~ 60 - 80 



S. Takami, K. Yabana and M. Matsuo* 
Graduate School of Science and Technology, Niigata University, Niigata 950-21 
* Yukawa Institute for Theoretical Physics, Kyoto University, Kyoto 606-01 



Abstract 



We study static non-axial octupole deformations in proton-rich Z = N nuclei, 
64 Ge, 68 Se, 72 Kr, 76 Sr, 80 Zr and 84 Mo, by using the Skyrme Hartree-Fock plus 
BCS calculation with no restrictions on the nuclear shape. The calculation 
predicts that the oblate ground state in 68 Se is extremely soft for the Y33 
triangular deformation, and that in 80 Zr the low-lying local minimum state 
coexisting with the prolate ground state has the Y32 tetrahedral deformation. 
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Spontaneous intrinsic deformation in the ground and excited states is one of the most fun- 
damental properties of nuclei. In addition to the well-established quadrupole deformation, 
octupole deformation violating reflection symmetries has been attracted much experimental 
and theoretical attentions. Recently, a pear-shape deformation, i.e., an axially symmetric 
Y30 octupole deformation superposed on a prolate quadrupole deformation has been estab- 
lished in light actinide (Ra-Th isotopes around A ~ 220) and lanthanide (Xe-Ba isotopes 
around A ~ 136) mass regions M. More exotic shapes such that violate both reflection 
and axial symmetries are known in light p and sci-shell nuclei, where the exotic shapes are 
caused by a cluster correlations 0, like as the equilateral triangular 3a configuration in 12 C. 
Besides cluster models assuming such configuration, mean-field models have also been 
discussed them in terms of non-axial octupole deformations. It is of great interest to see 
whether the non-axial exotic octupole deformations realize in heavier systems. Possibilities 
of the non-axial static octupole deformations have been argued theoretically for superde- 
formed nuclei |||7]] and for the light Ra region ||, for which however there seem to be no 
experimental support so far. 

Proton-rich nuclei in the A ~ 60 — 80 mass region are advantageous candidates for a 
search of the exotic reflection asymmetric deformations. For nuclei with neutron and proton 
numbers around 30-40, octupole correlation can be expected because of Aj = 3 coupling 
between the 2p 3 / 2 and lgg/2 orbitals in the major pfg-shell |5,10,lJ. Furthermore, the shell 
effect which is the important origin of the deformations is generally large in comparison with 
those in light lanthanide and actinide nuclei. Especially in proton-rich Z = N nuclei, both 
proton and neutron configurations coherently operate to develop a static octupole deforma- 
tion although there is no available experimental information such as E x (3i) and B(E3) in 
proton-rich Z = N nuclei around A ~ 80, except for 64 Ge [11]. The previous theoretical 



studies based on the shell correction method suggested that the octupole instability is not 
strong enough |f],n],E|| to develop a static octupole deformation in this mass region. How- 
ever, it should be noted that these models considered only the axially symmetric deformation 



PJTl|j, or limited combination of non-axial deformation parameters |fL2j| . 

In this letter, we will show that the reflection asymmetric shapes violating axial symmetry 
are more favored in Z = N nuclei in the mass A = 60 — 80 region than that with axial 
symmetry. We perform a fully self-consistent mean-field calculation for even-even nuclei 
64 Ge, 68 Se, 72 Kr, 76 Sr, 80 Zr and 84 Mo by means of the Skyrme Hartree-Fock plus BCS 
method with three-dimensional (3D) mesh representation |13|]. In contrast to the previous 
calculations with 3D mesh representation, we impose no requirement of symmetry on the 
solutions to allow arbitrary nuclear shapes. 

In the present method, quantities such as single-particle wave functions and densities are 
described on a 3D Cartesian mesh within spherical box. The radius of spherical box and 
the width of the 3D mesh are set to 12 fm and lfm, respectively. Imposing the constraints 
which diagonalize the mass inertia tensor, we make the principal axes coincide with x, y 
and z axes of the mesh. We adopt the Skyrme III interaction, which has been successful in 
describing systematically the ground state quadrupole deformation in proton and neutron 



rich Kr,Sr,Zr and Mo isotopes |13[ and in a wide area of nuclear chart ||14|| . It is also able 
to describe the axially symmetric octupole deformation in the Ra-Th region [15|| . As for 
the pairing strength of proton, we use the same parameterization G p = 16.5/ (11 + Z) MeV 
as in Ref. [O] together with the same truncation of the single-particle space. The neutron 
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pairing strength is taken the same as G p ||16|| . To characterize deformation of the obtained 
solutions, we have calculated the mass multipole moments, 

_ 47r($|S^rjX, m (»)|$) 
a/m = , (m= -I, ■■■,1), (0.1) 

where A is the number of nucleon and R = 1.2v4 1//3 fm. Here X[ m is a real basis of the 
spherical harmonics, 

Xio = Yio, 

Xl\m\ = —7={Xl-\m\ + ^-|m|)> 



-i 
V2' 



Xi-\ m \ - -j=(Xl\m\ ~ Yl* m \), (0-2) 



where the quantization axis is chosen as the largest and smallest principal inertia axes for 
prolate and oblate solutions, respectively. For the quadruple moment, we also use ordinary 
notation, i.e. , a 20 = /3cos7, V2a 2 2 = fisiwy, mapped in the (3 > 0,0 < 7 < n/3 
sections. To represent magnitude of the octupole deformation, we define 

/^(E^L) 1 , ^m = (a 2 3m + al_ m ^ (m = 0,l,2,3). (0.3) 

m=— 3 

For nuclei around A ~ 80, existence of three local energy minimum states showing oblate, 
nearly spherical and prolate solutions has been reported in the SHF+BCS calculations 
|13|]14T| . To search all minimum states close energetically to the ground state, we gener- 
ate initial states by solving a deformed Wood-Saxon potential model. The five initial states 
with different quadrupole deformations are used : (1) (3 = 0.7,7 = 60°, (2) (3 = 0.3,7 — 60°, 
(3) f3 = 0.0,7 = 0°, (4) (3 = 0.3,7 = 0°, (5) (3 = 0.7,7 = 0°. For all initial configurations, 
the octupole distortion, a 3m = 0.1 (m = —3, • • • , 3), is added. 

In table |, we summarize the quadrupole and octupole deformation parameters of the 
obtained solutions and the excitation energy of the local minimum states, where we do not 
list states higher than the fourth minimum. The solutions are classified into three groups, 
oblate, spherical and prolate, according to their quadrupole deformations. As is seen in 
table |, the calculated quadrupole deformation of the ground states varies significantly as 
changing the mass number. For Z = N = 32 ( 64 Ge), the calculation results in the triaxial 



deformation, which is in agreement to the experimental indications (TT|. It is also seen 
that the ground state quadrupole deformation changes suddenly from the moderate oblate 
ones in 68 Se and 72 Kr to the strong prolate ones in 76 Sr and 80 Zr. Sudden change of the 
quadrupole deformation qualitatively agrees with the observed systematics of the first excited 
energy levels (considered to be J n = 2 + state) from 64 Ge to 84 Mo [I7||. For all the calculated 



elements, we find the low-lying local minimum states with different quadrupole deformations, 
which is a characteristic of shape coexistence nuclei. 

The octupole deformations violating the axial symmetry are found in the ground state 
or local minimum states in all nuclei, except 64 Ge (where the obtained f3% = 0.01 is not 
sizable). Among them, the second minimum state of 80 Zr shows the largest octupole defor- 
mation of /?3 = /?32 = 0.24 without having a quadrupole deformation. The density profile 
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of this solution shown in Fig. [I] (a) indicates a tetrahedral deformation, which violates the 
both reflection and axial symmetries, but obeys the symmetry of the point group T^. Figure 
|2] (a) shows the potential energy surfaces of 80 Zr with respect to the a 30 , a 31 , a 32 and a 33 
deformations, which are calculated under the additional constraints of the mass octupole 
moments. The potential energy surface of the a 3 2 deformation has the minimum point at 
ct 3 2 = 0.24 which corresponds to the calculated lowest minimum, and the energy gain mea- 
sured from the spherical solution is as large as 0.71 MeV. This energy gain reduces the energy 
difference between the spherical and strongly prolate deformed solutions from 1.61 MeV to 
0.90 MeV. Octupole instability towards the a 3 2 direction (the tetrahedral deformation) is 
quite contrasting to the other types (a 3m (|m| ^ 2)) of the octupole deformations. 

Instability of the spherical configuration at Z = N = 40 for the tetrahedral deformation 
can be ascribed to the shell effect formed in the potential having the Td symmetry. In Fig.|2| 
(b), we display the neutron single particle energies as a function of the tetrahedral defor- 
mation parameter a 3 2- As developing the tetrahedral a 3 2 deformation of 80 Zr, the orbitals 
stemming from 2p 3 / 2 and 2pi/ 2 decrease in energy and those stemming from lgg/2 increase 
with holding high degeneracy of orbitals. The sub-shell gap at nucleon number 40 enhanced 
by addition of a 3 2 distortion field stabilizes the strongly tetrahedral deformed solution. It 
is known that high degeneracy of irreducible representation of the symmetry tends to 
produce a significant bunch in the single particle level spectrum as has been demonstrated 



for electrons in a metallic cluster potential by Hamamoto et. al. [|T^,|T^]. This tendency 
exists in nuclear potential with spin-orbit force. Appearance of the tetrahedral deformation 
due to similar shell effect will not be confined in this neutron/proton number as discussed 
by Li and Dudek for light actinide isotopes 0. It is also noted that the tetrahedral octupole 
deformation competes with the pairing correlation which favors spherical deformation. The 
energy gain caused by the tetrahedral deformation is found from calculation to be more 
than 2 MeV when we neglect the pairing. The tetrahedral minimum remains even with 
30% increase of the paring strength (G p ,G n ~ 0.4 MeV). It should be mentioned that the 
measured excitation energies of the fist 3~ levels in Ge and Se isotopes have the minimum 
points at N = 40 [^Ofl , which may be a fingerprint of octupole instability. 

The other types of non-axial octupole deformation are also found. The 68 Se is the most 
noticeable since it has the large octupole deformation (/? 3 = 0.15) in the ground state. 
As shown in the density distribution plotted in Fig. 1(b), it has V 33 triangular distortion 
superposed on the oblate quadrupole deformation, which obeys the D^h symmetry seen in 
the regular triangular prism shape. Although the minimum is not as deep as the tetrahedral 
deformation in 80 Zr, the potential energy surface is quite flat up to a 33 ~ 0.2 as shown in 
Fig. 3 (a). It should be noted that octupole instability emerges only for the a 33 direction. 
Instability toward the triangular deformation is systematically seen for the excited oblate 
minima in neighboring 72 Kr and 76 Sr. 

Instability of the oblate states toward the triangular Y" 33 deformation can also be related 
to the single-particle shell structure formed in the oblate deformed potential. Figure |] 
shows the neutron Nilsson diagram as a function of quadrupole deformation obtained in 
the constrained SHF+BCS method, in which axial and reflection symmetries are imposed. 
In the oblate configuration of 68 Se, the N, Z = 34 Fermi surfaces are located between the 
positive parity orbitals with Q = 9/2, 7/2, 1/2 stemming from the l<?g/2 and the negative 
parity orbitals with Q = 3/2, 1/2 arising from the 2p 3 / 2 (those just below the Fermi surface, 
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See Fig. 4). Among the possible couplings associated with the octupole deformations, the 
AQ = 3 coupling between the positive parity SI = 9/2 and negative parity SI = 3/2 orbitals, 
and also the one between the positive parity Q = 7/2 and negative parity Q = 1/2 orbitals 
have the smallest energy difference, and give enhanced softness toward the triangular Y 33 
deformation. 

In the prolate solutions of 68 Se and 72 Kr, the non-axial octupole deformations with finite 
value of /?3i (so-called 'banana' shape) are also found although the magnitude of the deforma- 
tion is small. The associated potential energy surface for the octupole deformations, shown 
in Fig||(b), indicates that the a 3 i direction is soft while the other types including a 30 pear- 
like deformation are not. It is noted that the positive parity orbital with Q = 3/2 (stemming 
from 1(79/2) and the negative parity Q = 1/2 (2p 3 / 2 l/i/2) located near the N, Z = 34 — 36 
Fermi surface (see Fig.4) are responsible for the Y31 deformation. Although there exists 
closely-lying AQ = pair of Q = 1/2 (lg 9 / 2 ) and SI = 1/2 (2p 3 / 2 l/i/2) orbitals, the matrix 
element of the AQ = octupole operator r 3 Y 30 between these orbitals is expected to be 
small, provided these orbitals have good asymptotic Nilsson [440] ~ and [310]| components, 
respectively. 

In summary, the self-consistent mean-field calculation predicts that 80 Zr has the low-lying 
coexisting state showing the spontaneous Y32 tetrahedral deformation and the oblate ground 
state in 68 Se is extremely soft with respect to the Y33 triangular deformation. Appearances 
of the non-axial octupole deformations are governed by the single-particle shell structures 
at N, Z ~ 40 and N, Z ~ 34 for the tetrahedral and triangular deformations, respectively. 

The non-axial octupole deformations discussed above will accompany specific spec- 
troscopic properties of excited states. The intrinsic tetrahedron shape having the Td 
symmetry leads to a series of levels with I 71 = + , 3~, 4 + , 6 + , 7~, • • •, which follow the 
E(I) — E(0) ~ /(/ + 1) relation and are connected with strongly enhanced E3 transi- 
tions H . In 80 Zr, for example, a low-lying excited J n = 3~ state is expected, and one 
could observe a strong E3 decay from a 3~ state to a low-lying + state associated with the 
tetrahedral deformed second minimum. The oblate deformed ground state of 68 Se which is 
extremely soft towards the Y33 direction will accompany a low-lying 3~ excited state with 
high collectivity. Furthermore, the intrinsic triangular shape having the D 3 h symmetry may 
form the characteristic band structure, that is, the ground band I + = + ,2 + ,4 + , ■ • ■ and a 
K w = 3" band I = 3", 4", 5~, • • • as seen in 12 C [§. 

The authors would like to thank Prof. K. Matsuyanagi for helpful discussions. Numerical 
computation in this work was mainly performed at the Yukawa Institute Computer Facility. 
Part of numerical calculations were also performed on the FACOM VPP-500 supercomputer 
in RIKEN. 
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FIGURES 



FIG. 1. Density distributions of proton in the xy, yz and zx planes where x, y and z axes 
represent the principal inertia axes, (a) and (b) show those of the second minimum state of 80 Zr 
and the ground state of 68 Se, respectively. 

FIG. 2. (a) The potential energy surfaces of 80 Zr with respect to the different types of the 
octupole deformations, where the energy is measured in relative to the spherical solution. The 
potential energy is calculated as a function of a^ m (m = 0,1,2,3) by imposing the constraints of 
[3 = 0,7 = 0° and a^ v = ( v ^ m ). (b) The single particle energy of neutron as a function of 
the tetrahedral 032 deformation. 

FIG. 3. Potential energy surfaces with respect to the different types of octupole deformations, 
calculated for the oblate ground state (a) and the second minimum prolate state (b) of 68 Se. The 
quadrupole deformations are set to (3 = 0.25,7 = 60° and (3 = 0.25,7 = 0° for the ground and 
second minimum states, respectively. 

FIG. 4. The neutron single-particle levels for 80 Zr as a function of the quadrupole deformation 
parameter calculated with the quadrupole constraint and the axial and reflection symmetries. 
For each orbitals, we put the value of f2, the projection of the angular momentum along the 
symmetry axis. The arrows indicate the Afl = 3 coupling associated with the triangular Y33 
deformation as discussed in the text. 
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TABLES 





Oblate 




Spherical 




Prolate 






g.s. 






0.62 


4.00 


64 Ge 


0,7 = 


0.27,25° (T) 




0,7 


= 0.24,6° 


0,7 = 0.38,0° 




03 = 


P33 = 0.01 




03 


= 0.00 


03 = 0.00 






g.s. 






0.32 


2.42 


68 Se 


0, 7 


= 0.25,60° 




0,7 


= 0.25,0° 


0,7 = 0.40,18° (T) 




03 = 


P33 = 0.15 




03 = 


031 = 0.06 


03 = 031 = 0.02 




g.s. 


1.12 






1.74 




72 Kr 


0, 7 = 0.34,60° 0,7 


= 0.27, 58° 




0,7 


= 0.42,1° 






03 = 0.00 3 = 


033 = 0.05 




03 = 


031 = 0.03 








2.58 


3.25 






g.s. 


76 Sr 


0, 7 


= 0.13,60° 


0,7 = 0.02,0° 






0,7 = 0.49,0° 




03 = 


033 = 0.16 


/?3 = 032 = 0.12 






03 = 0.00 






1.58 


0.90 






g.s. 


so Zr 


A 7 


= 0.20,59° 


0,7 = 0.00,0° 






0,7 = 0.50,0° 




03 = 


32 = 0.04 


03 = 032 = 0.24 






03 = 0.00 






g.s. 


0.24 






0.85 


84 Mo 


A 7 


= 0.20,56° 


0,7 = 0.05,60° 






0,7 = 0.64,0° 







j = 0.00 


03 = 030 = 0.13 






03 = 0.00 



TABLE I. The ground states and the local minimum states obtained in the present SHF+BCS 
calculation. The energy difference (MeV) between the ground state and the local minimum state 
(the ground state is denoted as g.s.), the quadrupole and octupole deformation parameters are 
shown, where label (T) stands for the triaxial deformation. For the solutions showing the octupole 
deformation, the most dominant component of 03 m is also presented. Each solutions are classified 
into the three groups, oblate, spherical and prolate by their quadrupole deformation parameter, 
except for the ground state of 64 Ge. The ground state of 64 Ge which shows the triaxial deformation 
is classified into the oblate group. 
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